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Abstract 

Nonlinear electrodynamics model in hypercomplex form is con- 
sidered. Its linearization around a solution is obtained. The appro- 
priate problem for linear waves around static dyon solution (SDS) 
of Born-Infeld electrodynamics is investigated. Two types of wave 
scattering on SDS are considered: dissipative (with momentum 
transmission from plane wave to SDS) and non-dissipative (for SDS 
imbedded to an equilibrium wave background). Resonance phe- 
nomenon in the problem is discovered and some resonance frequen- 
cies are obtained by using a numerical method. The form of reso- 
nance wave modes are discussed. The sum of a plane wave (as the 
elementary component of the wave background) with one resonance 
mode is considered. The appropriate energy density is investigated 
at infinity. The averaged energy density is demonstrated to have 
the term proportional to inverse radius. This fact allow to consider 
such field configurations as the cause of gravitational interaction, 
taking into account the effective Riemann space effect discovered in 
my previous works. A behavior of the linearized solution at origin of 
coordinates and the problem beyond the linearization are discussed. 
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1 Introduction 



Let us use the accordance between electromagnetic particles and space 
localized solutions of suitable electrodynamic field model. These model 
solutions will be called particle ones. 

Consider static particle solution of some electrodynamics model. Such 
a solution is the space-localized static field configuration satisfying model 
equations. For the usual linear electrodynamics we have a static particle 
solution in the form of Coulomb potential that corresponds to electrically 
charged particle. There are appropriate solutions for nonlinear electrody- 
namics. But in the case of nonlinear model the static solution will modify 
conditions of propagation for small additional electromagnetic waves. That 
is the space with a static solution looks like heterogeneous and anisotropic 
optical medium for these waves. The space-localized heterogeneity of a 
special kind, as with a static particle solution, must give a scattering ef- 
fect. This scattering may have a set of resonance frequencies determined 
by a specific size of the heterogeneity region. 

Dyon particle solution has point singularity with electric and magnetic 
charges. Here I consider a nonlinear model known as Born-Infeld electro- 
dynamics and its dyon static solution [3]. This approach can help us to 
advance in solution of some fundamental physical problems connected with 
particle physics. 

A similar problem for scattering of light by electrically charged parti- 
cle solution of Born-Infeld electrodynamics was considered by Schrodinger 
[H [10] , where the solution was seeking in the form of power series in fre- 
quency of light. But this form of the solution can be appropriate only for 
small frequencies. It should be noted that the corresponding mathematical 
problem is difficult for analytical investigation. 

The present state of the art in computational modeling allows to suc- 
ceed for the problems, which still remained unsolved just because of its 
mathematical difficulty. 

In the present work I use a method for investigation of the problem 
based on numerical calculations. This method is suitable for arbitrary 
characteristics of the electromagnetic waves. 
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2 Mathematical model 

for nonlinear electrodynamics 

Nonlinear electrodynamics in hyper complex or Clifford number form [1] is 
the most convenient to use. Here I consider the known representation for 
nonlinear electrodynamics, in which electric and magnetic inductions are 
used as dependent variables (see also my paper [3]). We have the following 
hypercomplex form for the basic equation without singularities [7J : 

d o Y + 0<->Y + #X\Z = O , (2.1) 

where <9 M = dj dx^ (Greek indices take the values 0, 1, 2, 3) and d = b b* <9j 
(Latin indices take the values 1, 2, 3) is the operator of space differentiation. 
Here b M are hypercomplex numbers appropriate to basis vectors. Symmet- 
rical <S> and asymmetrical A products are constructed from noncommuta- 
tive product for the hypercomplex numbers. For coordinate systems with 
the components of metric moo = — 1 and mo« = we have § = lb 4 di, where 
lb* are basis bivectors. At the present paper I use a coordinate system of 
this type. 

We have by definition 

Y = D + *B , Z = E + *H, (2.2) 

where i is hyperimaginary unit, D = D l fo>i and B = B l fl»i are electric 
and magnetic inductions, E = E l fl»j and H = H l Ifoj are electric and 
magnetic intensities. The quasi-bivectors of electromagnetic induction Y 
and intensity Z are invariant only for space transformations of coordinate 
system but not for transformations affecting time. The hyperimaginary 
unit i is permutable with do, Y, Z, and d. Thus for this case it can be 
used as conventional imaginary unit. The prefix "hyper" can be omitted 
in some words of the present paper. 

Symmetrical product is connected with scalar product for space vec- 
tors, for example D B = D • B, and asymmetrical one is connected with 
their vector product, for example DX\B=iDxB. But symmetrical and 
asymmetrical products have more wide sense for hypercomplex numbers. 
For the details see also my paper [7]. 

Equation ([2]) must be supplemented by the constitutive relations, con- 
necting electromagnetic inductions with intensities: D = D(E,B), H = 
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H(E,B). This connection can be expressed with the help of one scalar 
function £ = £(D, B), which is an energy density for the model. We have 
for the components 

d£ d£ 

In view of ( 12. 21) we can consider a function £(Y, *Y) instead of £(D, B). 
Here *Y designates hyperconjugate hypernumber for Y, realized with sub- 
stitution of * by — * in canonical representation ( 12. 2ft . 

We can consider strong or Frechet derivatives [8] of some function with 
respect to Y and *Y, where the bivector module |Y| = VD 2 + B 2 is con- 
sidered as the appropriate norm. These strong derivatives are linear opera- 
tors working on a function placed on the right of the operator. The strong 
derivative appears, in particular, in differentiation of function for vector 
argument. For example, 

d _ d£ dY d£ <9*Y 

■ £{Y, ) — — — : — — + — — : — — , (2.4) 



Qjyi v ' / dY ' dD l d*Y ' dD i 

where the symbol : is introduced to distinct the operation from the non- 
commutative product. If the strong derivative exists, then there is also a 
weak derivative equal to strong one. The weak derivative is defined by the 
regular way: 



Sf < Y > : Q = 



|fCY + eQ) 



OY 

According to (TOD, and (T23D we have 



(2.5) 

£ = 



8£ ~ 1 ~ d£ ~ 1 

^-:Y=— *Z0Y , ^-:*Y=—Z®*Y. (2.6) 
dY 8tt ' d*Y 8tt V ; 

For Born-Infeld model (see also [6]) we have the following expression 
for energy density (see, for example, my papei0 [3]): 



1 In the present paper the designations are somewhat changed. 
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where 

€ = ^l + 87f X 2 £ + (4tt) 2 x 4 V 2 , (2.7b) 

€ = — YO*Y = — (D 2 + B 2 ) , (2.7c) 
8tt 8ir 

?4 Y ^ = -7-( D » B )4( DxB ) > (2-7d) 

87T 47T 47T 

X is a dimensional constant of nonlinearity for the model. Here £ is energy 
density for linear electrodynamics, that is for |Y| — > 0. "P is momentum 
density both for linear and nonlinear electrodynamics. It is obvious 

POY = . (2.7e) 

Having a form for the energy density £(Y, Y), we can write an appro- 
priate nonlinear constitutive relation according to (12.31) and (12.2ft : 

Z = Z(Y, Y) . (2.8) 

The constitutive relation specifies a concrete nonlinear electrodynamics 
model. 

According to (12. 7p . (12. 3p . (12. 1\ we have the following constitutive rela- 
tion for Born-Infeld model in hypercomplex form: 

1 + 47T Y 2 V 

Z = — * Y . (2.9) 
€ 

This constitutive relation differs from one for the case where electric in- 
tensity E and magnetic induction B are considered as dependent variables 

HI- 



3 Linearization around some solution 

o 

Consider a solution Y for nonlinear electrodynamics model and a lineariza- 
tion Y around Y. Thus I find a solution to equation ( 12. lft in the form 

Y = Y + Y , (3.1) 

where Y is small: 

\ X Y\ «1 . (3.2) 
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Here x is considered as the dimensional constant for general electromag- 
netic field model with arbitrary constitutive relation (12.81) . 

Substitution of (13.11) into constitutive relations (12.81) and consideration 
(13.21) give the following linearization for the general case: 



Z = Z + — 
dY 

where Z = Z(Y). 

For relation (12. 9p we have 



: Y 



Y=Y 



dZ 



d*Y 



. *y" 



Y=Y 



(3.3) 



z = z + z , 

Z = E- 1 { Y + ( X 2 /2) 



YA*YAY+YA*YAY+YA*YAY 



• (3-4) 



- Z(^ZO*Y + *Z<»Y 

Here we can apply the following useful formula: 

(BAC) A A = B(AOC) -C(AOB) . (3.5) 



4 Static dyon solution (SDS) 

The simplest spherically symmetric static solution of Born-Infeld electro- 
dynamics has one singular point with electric and magnetic charges and 
called dyon [3]. This static dyon solution (SDS) can be written in the 
following hypercomplex form: 



Y 



G 



oo 



(4.1) 



where G = (§ + i & is electromagnetic charge, G is electric charge, G is 
magnetic one, are angle spherical bivector functions introduced in my 
paper |?J , ©{j = lb r is radial basis bivector, f = \f\x~ G| is a characteristic 
radius of the static dyon solution. 

G 

Introduce also a normed charge G = — -, which characterizes a ratio 



G 



between electric and magnetic charges. 
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5 Linearization around SDS 
and appropriate solutions 

Let us consider a time harmonic field configuration expanded in terms of 
the angle spherical bivector functions <§^ m [TJ: 

1 oo / 

Y = Gj2J2J2^r (+fYf e-^ z ° + _ f Yj?V wai0 ) , (5.1) 

j=-l l=\j\ m=-l 

where ©J™ = <p),_ ±f Yj m = ±$ rl ™(r), cu > 0. Here the multiplication 

by the normed charge G is used for a simplification of resulting formulas. 
The analogous representation is considered for quasi-bivector Z. 



5.1 Equations 

Substituting (15 .ip into (13.4)) . taking into consideration (14.11) . and using the 
dependent variables introduced in [7], I obtain 



flm 



rjlra 
±1^0 



±2 

where 



1/7 1 



Vr 4 + f 4 
Vr 4 + f 4 



Vr 4 + f 4 



±f 1 + 



2r 4 



±f 1 + ~+~ 



-l) m T ^ 



L—m 



±f* 



2 (r 4 + f 4 ) 



if- 1 



T± X 



if 1 - 



-4 

2T 4 



±f 



Y lm - (-1) 



Tf 1 " 



if 1 -1 ' 



if 1 + — if 1 1 "T ±f I _1 , ±fl _ — ±fl 1 

rvim _l_ ryim I rjlm rylm _l_ rvim vim 

±f£ + — ±f ^i + ±1^-1 > ±f^- — ±f^l _ • 

According to [7] there is the following system of equations: 
(r 2 ±f Y^). r + i r + ±f Yf = , 



r 



lm 


im 



±f 



E Y< m ± (Z + l) ± £!? 








(5.2a) 

}• 

(5.2b) 
(5.2c) 

(5.3) 



(5.4a) 
(5.4b) 
(5.4c) 
(5.4d) 
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where (...) ;r is derivative with respect to radius r. 

To simplify relations (15.21) introduce again the following new variables: 



•\rlra _i_ -\/-lm 
±|f X p — +f 



£ p 31 V - 1 / f 1 p 



ryim _i_ r^Zm if i \m *rvL—m 
±\1 L P — +1 L p ± \~ L ) -f^p > 

(5.5) 



where p = — , 0, +. 

Inverting relations (15.31) and (15. 5p we have 



+f 1 



\rlm 
-f 1 



•\rlrn 
+f* ±1 



-f 1 ±1 



"V~Zm _i_ ~\rlm 

2 V+|f* o "I" -if 1 o 



" 1)m r * y 



L—m 



* -Crl,—m 
-If 1 



rZ,— m 



+ |f J 



■ V"' 



"|f J 



=t 1 +|f* - 



L—m 



Using (15. 5p I obtain the following relations from (15. 2p : 



(5.6) 



r Zm 



Vr 4 + f 4 



If* + ' 



r^Zm ' 

+ 1^0 — / 4 _|_ ^4)3/2 



-£/-Zm 
+ |f Y 



rjlm 



-\rlm 

If - ' 



'Zm 


r 2 


+ - 


Vr 4 + f 4 


'Zm 




o - 


Vr 4 + f 4 


'Z»7i 


Vr 4 + f 4 






(13751): 





-£/-Zm 
-|f Y + ' 



-If 1 ) 



-I* 



Y 



Irn 



(5.7a) 
(5.7b) 
(5.7c) 



*rv^7TTl) T | f Y 



Zm 











r u; 



(r ±| ^_ m ) ; ,-ir^ ±|f V 
£ Y^ m + V/ (Z + l) ,^ = 



±|f 



(5.8a) 
(5.8b) 
(5.8c) 
(5.8d) 



Differentiating (I5.8al) and making necessary substitutions from (15.71) 
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and (15. 8p I obtain the following two equations for the functions ± | f YQ m (r): 

(d*_ 2(r 4 + 2f 4 ) d 
\ dr 2 r (r 4 + f 4 ) dr 

(5.9a) 

}(r Hf *S")=0. 

(5.9b) 

These equations can be written also in the form 

r 2 +|f Y^) =0, (5.10a) 

r 2 Hf Y[, m ) = . (5.10b) 

We can introduce the following dimensionless radius and frequency: 
r 

r = - , Cj = ojr , ijjr = ujr . (5-11) 
r 

The appropriate "dimensionless" equations are obtained with the substi- 
tution f — > 1, uj —> uj, r — > r into equations (15.91) and (I5.10p . 

The operators of equations (I5.10p coincide with operators of equations 
considered by Schrodinger [9l [TOj . But dependent variables in Schrodinger's 
approach was electromagnetic potential components. 

If we have the component ± |^o m fr° m solution of equation (15. 9p or 
(I5.10p . then we obtain immediately all components ± ^ l ± and ±\ft™ with 



|(r +|f Y^)=0, 



+ 



2H 



dr 2 (r 4 + r 4 ) dr 



d 

-r- + 



2 Z(Z + l)r 4 + 2f 4 

j.2 ^4 _|_ 



d^_ 

dr 2 
dr 2 



2f 



r (r 4 + r 4 ) dr 
2f 4 



+ 



' 2 Z(/ + l)r 



21 



r (r 4 + r 4 ) dr 



d 

-r- + 



a; 



(r 4 + f 4 ) 
/(/ + l)r 2 
(r 4 + f 4 ) 
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the help of (J£2) and (Oil : 

+ | f I — <y + C, + , _| f I o — W_ 



■^yim 2 jL. '\rlrn 

l f + ~ /777~TTT^ ' "l f 



r 



I A) " (r ~4 + 1)3 /2 > -|A> - v ^4 TI 

/ ,/, , ,x -o ' -|A 



?i m « (f 2 cr ) „ ©^i? (r 2 cr).. 



(5-12) 

where © l ± are complex constants, e± J = C± w (^)- 

5.2 Solutions in asymptotic power series form 

The coefficients of equations (15.101) contain the characteristic function (1 + f 4 ) -1 
which can be expanded in asymptotic power series near zero and near in- 
finity: 

oo oo 

Y^-i = E(- 1 ) nf4n = E(- 1 ) n+1 ^ 4n • ( 5 - 13 ) 

n=0 n=l 

The convergence domain for the first series in ( 15.131) is f < 1, and for 
the second one it is f > 1. Thus we can not expect a convergence of 
appropriate power series representation for the solutions for all f values. 
Let us consider two asymptotic expansions. 

At first consider solutions to equations (15.101) in a generalized power 
series form near point of origin: 



(> Cf ) = f ( 1 + f; IB? r n J , (5. 14) 



n=l 



where C± w (^) are functions introduced in (15.121) . p and [n ^ 1) must 
be obtained. 
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By direct substitution of (I5.14p into (15.101) we have p = — 1 or p = for 
(I5.10al) and p = or p = 3 for (15.10bj) . But the cases p = — 1 for (I5.10al) 



and p = for (15.10bl) give rise to the singularity of types + | f Y 



Irn 



and 



+|f J 



r lm 



that give the appropriate divergence of an energy integral 
(see also (12.71) ). Thus we must take p = for (I5.10al) and p = 3 for 
(IBTObj) . The coefficients in flSTTij) can be obtained easily with the 
help of recurrence formulas given by equations (15.101) with the substitution 
of (I5.14p . Let us write the appropriate solutions containing several terms: 



Co 2 



6 



cD 4 



— 1(1 + 1) + 

20 v ' 120 



— r 



Cu 6 



I (I + 1) H 1 

2520 v ; 63 5040 



r^0 



) I ~2 ^ 

r {1 — r — 
10 



f 4 



— 1(1 
28 v 



' 14 + 280 



— r 



19cD 2 ,., _ 17cD 2 



u 6 



7560 



756 



15120 



(5.15a) 



(5.15b) 



where O^Z+o denotes a quantity having the same infinitesimal order that r n 
as r — > 0. 

Eqs. ( 15. 91) are converted to equation for radial spherical functions 
(introduced in [7j) when we take f = or r — > oo. These functions are 
connected with spherical Bessel functions. We have the following behaviour 
of the radial spherical functions at infinity r — > oo (see (8.29), (8.31) in 

my- 

(5.16) 



Thus we can seek an appropriate asymptotic solution in the following form: 



~2 A,u 



(5.17) 



n=0 



where ^q 1 * 1 7^ 0, the complex coefficients for n > must be obtained, 
9ft is hyperreal part. Here I take the real part of the series to have the 
concordance with series (15.151) (see (I5.27P below in subsection 15. 3j) . The 
coefficients ® l ^ L (n > 0) can be obtained easily with the help of recurrence 
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formulas given by equations (15.101) with the substitution of (15.171) . Let us 
write the solution with several terms in the following form: 

(5.18) 

where €l l ± = 0q of ( 15.171) . Here the first term within the square brackets 
is characteristic for the spherical Bessel functions. 

Series (15.141) . (15 . 151) must be considered as asymptotic one for suf- 
ficiently small value of f. Series (15.171) . (15.181) must be considered as 
asymptotic one for sufficiently great value of r. Connection between these 
two asymptotic expansions is characterized by the complex constants @ l ±. 
These constants must be obtained with the help of another methods. I use 
a numerical calculation for this purpose. 



i 



('±|]L ( _ 2lwr) -.' ± _L + 



5.3 Numerical calculation from infinity 

Let us seek one of two solutions for each of two equations (15.101) in the 
form 

-^y l ± {f) , i4(oo) = l . (5.19) 

Here the functions y l ±{f) have zonal index /, which can be omitted in the 
following expressions. 

Change the independent variable f to s = r~ l and introduce new un- 
known function 

z±{s)=y±{s- 1 ) ■ (5.20) 

Substitution of (15. 19j) with (I5.20p into (15.101) gives the following equa- 
tions: 



. / 1 tu\ , I (I + 1) - 2ts 3 u 



/ +2 f L±l£. ' !»); Jl±pl^ z _ = . (5.21b) 

V s (! + ^ 4 ) s 2 J - s 2 (1 + s 4 ) V ; 

These equations are investigated with the condition z±(0) = 1. 

The numerical calculations of (I5.2ip with defined value of / starts from 
the point s — si, where z±(si) is calculated with the help of the power 
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series at infinity (I5.18p . The used value is s% = -f- 10~ 8 and obtained 
results does not depend on this value. The calculations is stopped at the 
point s = S2- To obtain a proper exponent of s at infinity or f at zero for 
the solution, I take the value of S2 to be sufficiently large: S2 = 10 9 . The 
obtained complex numerical solutions are represented in the form 

z±(s) = \z±\exp[i(z±}} , y±(r) = \y ± \exp[i(y ± }} , (5.22) 

where \z±\ is absolute value and (z±) is argument of z±. 

The numerical calculation gives \z + \ ~ s and \zJ\ ~ 1 as s — > oo. Using 
also the representation of y ± as power series in f , obtained from equations 
(I5.10p . we have 



a 1 

^ = _± e ^+ + ... , yi. = aLe^-(l-ta;r) + ... , (5.23) 

where a l + = lim [f |yf|_(f)|], a l _ 4= lim [||/L(f)|], (fi± = lim [(^(f))]. 

Second linearly independent solution for each equation (I5.10p is ob- 
tained by complex conjugation of first one (I5.19p . Thus we have the fol- 
lowing general solution to equations (15.101) : 



V C± j = 2 



(5.24) 



where l ± and '& l ± are arbitrary complex constants. 

We can obtain directly the following. In general case expression (15.241) 
with (15.231) is appropriate to the forbidden exponents in (15. 14[) : p = — 1 
for the "plus" component and p = for the "minus" one. To have the 
admissible exponents p = and p = 3 accordingly, we must put 



'& l ± 



-l) l €J l ± e 2i ^ . (5.25) 
Comparing also (I5.24p . (15.251) with (15.151) I obtain 



a + on a_uj 6 



a l + uj a l _uj 3 ~' 2 



(5.26a) 

I) = - - . (5.26b) 
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In view of (I5.26al) relation (I5.25P gives '0± = *0 l ± . Thus we can write 
the following, instead of (15.241) : 

(?<g)=*{&±%^ ■ (5-27) 

It should be particularly emphasized that here we have the necessity of 
reality for the functions Because this I take the real parts in above- 
stated expressions (15.1 7ft . (15.181) . 

At infinity (r — > oo) Eqs. (15.91) take the form of equation for spherical 
Bessel functions. Let us represent the solution at infinity with the help 
of radial spherical functions $~. Comparing (I5.27P with the appropriate 
definitions for S~ [7] I obtain the following expression for r — > oo: 

Cdf = ^4 [cos(^L) - sin(^ ± ) 3Sy , (5.28) 

where \@ l ±\ and (@ l ±) are defined in ( I5.26bl) . S~ = % l {oo r), S\ur) are 
spherical Bessel functions, 5?$' is the spherical Bessel function of the first 
kind, SsS* is the spherical Bessel function of the second kind. Here, of 
course, u l denotes the Z-th power of Co. 

The amplitudes a l ± and phases <j) l ± for y l ±(0) (15.231) obtained numerically 
for the values u = 0.05 -f- 600 and I = 1 4- 10 are shown in Fig. [TJ These 
calculations give also 

for cu^O , (5.29a) 



for u^oo , (5.29b) 

where the coefficients a l ± in (I5.29al) are rapidly increasing with increase 
of the zonal index /. For example a 1 ^ m 0.7, 5l ~ 1.3, a+ ~ 9.1 • 10 3 , 
a 5 _ » 3.6 • 10 3 , a™ w 5.1 ■ 10 10 , a 10 w 1.3 • 10 10 . 



j _ ~ a ± 



to 



a_ = — 



4 



7T 

~2 
(/ + 1)tt 

2 
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Figure 1: The frequency dependence of amplitudes a ± and phases 4> ± for 
//' (0) (E23D. 
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6 Linear waves around SDS 

Any solution Y can be represented in the form Y = Y_i + Yo + Yi, 
where each component Yj includes the appropriate radial functions -t f Y^ m 
of expansion (15. ip . According to (15.121) with (15. 6p the components Y±i 
are uniquely determined by the component Yo. Thus we can analyze this 
component only. Let us write its form at infinity according to (15. ip . (15. 6p . 
(I5~T21) . and ([52HD: 



Y = G^ 



l in 




1=0 m=-l 



-\7-lm —t uj x c 
-f 1 e 



i -Crlm xujx 
"I" _ f I o e 



where 

+f 1 



-f 1 



1 



— [(3*" \8 l + \ (cos(^) m 1 ^ - sin(0 l + ) SsS^) 



2Q l f 



+ 



Vim I ypfi 



2u l r L 



(cos(& l _) - sm(8 l _) 9f«y] 
(cos{0 l + ) m 1 ^ - sm(8 l + ) 3$^) 
- W^ m (cos(#!_) - sin(@ l _) 3$^) 



(6.1) 



(6.2a) 



(6.2b) 



Here @ l ± are defined by (I5.26a[) . (g! 1 ™ are arbitrary complex constants. 
The constants must be obtained from additional conditions. 



6.1 Linear wave scattering on SDS 

The known approach to the problem for wave scattering implies that an 
incident or acting wave is given for all space, and a scatterer gives birth 
to a scattered wave. This approach, when applied to the problem under 
consideration, suggests that the field Y at infinity includes two part: 

Y = Y + Y , (6.3) 

where Y is a field of the acting (incident) wave and Y is an additional 
field (scattered or response wave), such that the field Y is a solution for 
problem (15.81) . (15.71) at infinity. 

The condition for distinguishing of incident and scattered waves can be 
their behaviour at infinity: 

|Y| = O].^ , |Y| = CC, . (6.4) 
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But it should be borne in mind that decomposition (16.31) with condition 
(16.41) is ambiguous. Really, this condition remain valid if a term of higher 

order than 1 (o 1 ^^) is subtracted from Y and added to Y. 

Consider the acting wave mode in the form of plane wave propagating 
on positive direction of x 3 axis with arbitrary polarization. This wave can 
be represented as the linear combination of two waves with mutually an- 
tithetic circular polarizations. Using designations introduced in my article 
[7J we have the following representations for this wave (see Eqs. (7.28) in 

my- 

Y = C® 1 1 \® + e tuj{rcos& - x ° ) + ®-e- luj{rcos *- x0) ] , (6.5) 

where u > 0, CD™ are angle cylindrical functions, relation x 3 = r cos$ is 
used, the complex coefficient is used for the concordance with (15. ip . Two 
complex numbers ® + and ©_ (with given 0) define a polarization ellipse 
and an amplitude of the wave. 

According to [7J (see Eqs. (8.46) in [7j) we have the following expansion 



(K{e 



1 „iur cos # 



oo 

i=i 



where 



/ - V y/Tif+l) (21 + 1) 

will be called spherical coefficients of the plane wave, 



r 



Y l - - 
Jf±i - 2 



l + l 



iuj 



i-i 



r / 



(6.6a) 

(6.6b) 
(6.6c) 

(6.6d) 



Using representation (15. ip for the field Y ( I6.5P and comparing with 
( 16.61) . we have the following coefficients for the case of a plane wave with 
elliptic polarization and propagating in positive direction of a; 3 axis: 



±f* j 



-1-2 



- 3 



(6.7) 



r—*r, w— >±a 
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where §™ is the Kronecker symbol: §™ = 1 for m = n, §™ = for m 7^ n. 
Here I use dimensionless coordinates and frequency (15. lip . ux° = ujx°, 
x° = x°/f. 

The field of plane wave expanded in terms of spherical functions ( 16. 5p . 
(16.61) includes divergent and converging elementary spherical waves [TJ: 

± @ l ™e-™ x ° , (6.8) 

where the sign index + or — corresponds to divergent or convergent wave 
accordingly. 

The representation for solution as the sum of incident and scattered 
waves is usual for problems of wave scattering on bodies (for electromag- 
netic waves see, for example, [HUT]). In this case a progressing plane wave 
is taken as the incident wave, and a sum of divergent spherical waves is 
taken as scattered wave. Thus in this approach the scattered wave contains 
the divergent spherical waves only. 

According to this approach we can write: 

Jt? = - g %g M + Rf ay , ( 6 .9) 

where the first term in the brackets corresponds to the plane wave according 
to (16. 7p . The second term corresponds to the divergent spherical waves 
(see (8.38) and below in [7J), where complex constants R l ™ considered 
as "reflection coefficients" of convergent spherical waves contained in the 
plane wave. 

Obtain radial functions ± | f Yo ±:L from (16.91) with the help of (15.51) . then 
compare it with the radial functions from (15.121) and ( 15.281) . As result we 
have 

<$£ u 2 \& l ± \ [cos(^) m 1 ^ - sin(0 l ± ) 3Sy = -&+ t («S^ + R^ 1 , 

(6.10a) 

©^u 2 \@ l ± \ [cos(@ l ± ) - sin(& l ± ) 3Sy = T*®- @i («S^ + *R l ±~ 1 Sj-) 

(6.10b) 

where % l ±a = 3?S^ + tOfS?^ = ±3?S^ + and used ^ = (-1)' 

the constants @ l ± are defined by (15.261) . 
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Because of linear independence between 9ft§^ and we have from 
(jfi.lOp a system of linear algebraical equations for (5±\ R l +, ©jT 1 , *R±~ 1 . 
Its solution is 

= * R i r i = t exp ^ sin ^ . ( 6 .n) 

For m 7^ ±1 the coefficients @± m = 0. The solution ± f YQ m at infinity for 
linearized equation system ( 15. 4ft is given by formula ( I6.2p with coefficients 
(5± m obtained here. 

The obtained solution is based on the intuitive assumption for absence 
of convergent waves in scattered wave mode. But mathematically it is not 
quite clear why we must exclude the convergent waves from the scattered 
wave, because the divergent and convergent spherical waves are complete 
set of functions for the representation of solution. However the approach 
considered above, when applied to the usual scattering of light on bodies, 
allow to obtain the observable effect of momentum transmission from light 
to the scattering body. 

This scattering can be called a dissipative one, because there is the 
permanent energy-momentum transfer from the incident plane wave to the 
divergent scattered one. 



6.2 Equilibrium wave background around SDS 

We can also consider an experimental configuration with an equilibrium in- 
terchange between acting waves and response ones by the scattering space 
region. This situation can takes place with SDS imbedded to a wave back- 
ground. In this case the appropriate thermalization process must be ac- 
counted for by a weak wave interaction caused by the nonlinearity of the 
model. This wave interaction is not considered here. 

The weak wave background can be represented as expansion in terms of 
plane waves with various directions of propagation and polarizations. Thus 
in this case we can consider also a response of SDS to one plane wave (16. 5p . 
but taking into account the equilibrium between acting and response waves. 
In this case the response mode Y in (I6.3P must include divergent and 
convergent spherical waves, such that there is no the permanent momentum 
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flux through closed surface containing SDS. This scattering can be called 
a non-dissipative one. 

The appropriate solution at infinity for one acting plane wave, consid- 
ered as a part of the weak wave background, has the following form: 



ygn = _ g M (±1)' ^ + ^ ^) ( (612) 



±f 



where i?^™ are some coefficients. 



As we see, expression (16.121) differs from ( 16. 9ft with change i?± S ±(I , — > 
R l ™ 3$^,. Thus we can consider relations (16.101) with this change. Here I 

use designations with tilde for the coefficients Solving the appropriate 
system of linear algebraical equations and using (I5.26bl) I obtain 

^ = -©+0,3'+ , @ n = -^® + l ^ l _ , 

= -*®_ ( §* + , = - *&_ l Q l _ , = R l r X = - tan(^) 

(6.13a) 

where Q l ± will be called the response functions: 



cos{® l ±) 



(6.13b) 



As we can see in (16.131) . the coefficients and M^ 1 for solution 

(I6.12j) have infinite values for cos(& l ±) = 0. The condition cos(& l ±) 



0, with ( I5.26bl) and calculations presented on Fig. Q], defines resonance 



frequencies for the case of non-dissipative scattering. This resonance is 
connected with the radius f of SDS. 

It should be noted that this resonance appears also in the case of dis- 
sipative scattering, described in subsection 16.11 where cos(& l ±) = is the 
condition for maximum of the reflection coefficients R l ™. In this case the 
coefficients of the solution are limited because of the dissipation effect. 

The numerical calculation of the response functions (§ l ± for I — 1, . . . , 10 
are shown on Fig. [2j The pointed peaks on Fig. [2] correspond to the 
resonance frequencies. 



20 




Figure 2: The frequency dependence for squares of the response functions 
(16.13bj) . The relation between the values of zonal index I and the type of 



dashed curve is the same that on Fig. [TJ 

The values of the resonance frequencies (bf n are presented in Tab. [TJ 
where n is a number of the frequency with fixed zonal index I. We can also 
range the resonance frequencies in ascending values. The ranged resonance 
frequencies up to Co J, x ~ 5.45 are presented in Tab. [2j As we can see on 
Tab. [1], to obtain the correct range for u ln > uf 1 we must have the 
calculation for I > 10. 

The numerical calculations give also the following asymptotic for uj — > 

oo: 

^ + ^(-l)'#V , , (6.14) 

where can be represented by the approximate formula 

-i 1.86 1.91 2.45 
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But if we exclude the value for I — 1 then the good approximation is 

#+«^ . (6.16) 

For I = 1 this formula give a less value then one obtained with the numerical 
calculation 2.40. 



/ 




"in 


1 






2 






3 


2.45 


4.32 


4 


3.53 


6.97 


5 


3.52, 5.48 


4.51, 10.45 


6 


4.07, 7.56 


5.55, 14.66 


7 


4.42, 5.79, 10.04 


5.48, 7.20, 19.59 


8 


4.80, 7.06, 12.87 


6.16, 9.03, 25.23 


9 


5.12, 6.80, 8.62, 16.05 


6.48, 7.50, 11.12, 31.59 


10 


5.45, 7.73, 10.38, 19.58 


6.95, 8.86, 13.45, 38.65 



Table 1: Resonance frequencies 



Q 




/ 


n 


± 




1 


2.45 


3 


1 


+ 




2 


3.52 


5 


1 


+ 


1.44 


3 


3.53 


4 


1 


+ 


1.00 


4 


4.07 


6 


1 


+ 


1.15 


5 


4.32 


3 


1 




1.06 


6 


4.42 


7 


1 


+ 


1.02 


7 


4.51 


5 


1 




1.02 


8 


4.80 


8 


1 


+ 


1.07 


9 


5.12 


9 


1 


+ 


1.07 


10 


5.45 


10 


1 


+ 


1.06 



Table 2: Ranged resonance frequencies 
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We can consider a special case of the wave background with zero mo- 
mentum density. In this case the wave background can be represented as 
the sum of standing plane waves with various directions of propagation and 
polarizations. 

In this connection consider the linear response of SDS to a standing 
plane wave on x 3 axis with arbitrary polarization. This wave has the 
following form: 



Y = C 



+ (®\e~ lujx3 + ®z\e tujx ^e lujx ° , (6.17) 

where u > 0, the coefficients ® + and (P_ (with given G) define a polariza- 
tion ellipse and an amplitude of the standing wave. 

Taking the conjugation for f)6.6al) and taking into account the relations 

= (-1)' l and m l ™ = (-l) m I obtain 

oo 

(0; = i e -»«r cob* = _^_ 1 y+i 0^-1-2^1,-1 _ ( 61g ) 
1=1 

Thus for the case of standing wave (I6.17P we have the following coeffi- 
cients of representation (15. ip instead of expressions (16 .7p : 

Jtf = -[^ + ex(-l) m ]^0ia;- M ^}|^ iW _^ . (6.19) 

The coefficients of the solution for this case are 

@+ = -[&+- "©-] @i f + , @ l l = ~ [©+ + (-1)' *&-] ©i C , 

o 

(S^ 1 = [(-1)' &+ - "©_] 0i Q l + , 1 = - [(-1)' &+ + "©_] 0/ f [_ • 

(6.20) 



6.3 Resonance wave modes 

The effect of wave resonance connected with two-point boundary problem is 
well known. In this way the boundary condition on two resonator walls (at 
two points of appropriate normal coordinate) give rise to discrete frequency 
spectrum for waves within resonator, that is between the walls. 
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But some two-point boundary conditions can be noticed also in the 
problem of the scattering on SDS. In this case we have the boundary con- 
ditions at the singular point and at infinity The convergence of energy 
integral at the origin of coordinates is the first boundary condition. The 
condition at infinity is that the resonance mode has a form which is not 
representable by regular harmonic waves of linear electrodynamics. 

The last sentence implies that the resonance wave modes at infinity 
include only imaginary part of the radial spherical functions SS^, but 
regular spherical waves in linear electrodynamics include the appropriate 
real part only. 

At the resonance frequency uj = ujf n we have cos(®±) = 1, sin(®±) = 
in (I5.28P and the resonance wave modes have the following form: 

1 I x-W I 

S lnm = Yl if^^e™^ , where ± ^ l n -> ' * 1 _ SSj.± as r - 

(6.21) 

The resonance wave modes at infinity are orthogonal to the regular 
acting wave Y in the following sense: 

ri+27r/u) 7t 2tt 

Jim Jf 2 dr j&d Jdip (±Y W <J> ^ =0 , (6.22) 

fi oo 

This orthogonality is connected with the evident relation 

.Km J KSiS$if 2 df= J U[%- l - x e™*] ^[%- l ~ x B %ar \ dr = , 

(6.23) 

where expression (1 5 . 1 6 [) is used. 

By using (I5.16P the boundary condition for the resonance modes at 
infinity can be written also in the form: 

Jginm ^ 9«[ t -'-iexp(tw^r)] as f -> oo . (6.24) 



OO 
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7 Wave solution at infinity 
and gravitational interaction 

There is a possibility for unification of electromagnetic and gravitational 
interactions in the framework of the nonlinear electrodynamics because of 
an effective Riemannian space effect [21 El E] ■ This effect implies that the 
energy-momentum density tensor of electromagnetic field creates the ef- 
fective Riemannian space for propagation of weak quick-oscillating waves. 
This effect may explain the gravitation, if an expression of averaged en- 
ergy density for electromagnetic particle contains C^oo t erm ( see 0> the 
averaging is needed for cancel quick-oscillating terms). But for the static 
solution (14.11) we have only 0£-+oo behavior. Here I investigate this prob- 
lem for representation (16. 3p in the case of non-dissipative scattering for 
equilibrium wave background, stated in subsection 16.21 

According to (I5.28P and (I5.12p we have the following expressions for 

the components of Y at r > 1: 

±|f " ^ ~ ' " + 1) f ' ±|fY " " " y/l(l + 1) 

(7.1a) 

where 

Using (JESD, (17151) . fl5TT6|) . and (8.34a in [7J) I obtain 
Hi a, 

= J^^y [Tfi9f(*-'- 1 e wr ) -»S(r ! e- r )] + , 

-m , *&l,-m 

= (-l) m ^± ZJ ^=- \±Co S e™ r ) + 1 S ft" 1 e^ r )l 

1 ±X V 7 4f^ (J + l) 



-f 



+ (Coo > ±f Y r = <Coo • (7-2) 

Here I consider the case when the value of frequency Q is near to one 
of its resonance values: Cj ~ Q q (see Table |2]). In this case tan(0 l ± ) — ► ±00 

and ^f^ 1 — > ±00 for one value of I = l q and definite sign index + or — . 

Thus we can consider in a comparative sense that ~ for I 7^ Z 9 and 
the appropriate sign index + or — . 
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At infinity the energy-momentum density tensor has the form which 
is appropriate to linear electrodynamics. In particular, the energy density 
has the form (I2.7cl) . Thus I consider expression & for Y ( 16. 3p , where I seek 
the terms having C^->oo behavior for an averaging in the space-time period 
2n/u. 

The four-dimensional averaging for a function Q(x^) near some point 
x^ can be defined 

{Q)l/z 4= 77^4 / « dtf , (7.3) 
{2-k/uj) J 

where W is the four- volume, such that x M — tc/u ^ x'^ ^ x M + tc/u. Here 
we have the dependence {Q)L/u, = {Q)L/u, 

Substitution of (16.31) for (12.7cD allows to write the following expression 
for the averaging energy density of the solution at infinity (16. 3ft : 

PJir/a = W + * , V) + aCoo , (7.4a) 

where 

W=-1(y«>yT , (7.4b) 
cross terms with Y and Y give 

(y<E>¥+ ! Y<=>^) = I^(f,^, ¥ ,) + O^ 2 oo . (7.4c) 
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Here $^(f, tf, </?) = O 1 ^ is a function which must be obtained for each I 
and each sign index. It is evident that $i (r, ir/2,<p) = because the plane 
wave do not vary in the flat x 3 = and the averaging of (15.161) gives zero 
for C^oo terms. Here I obtain values $r (f, 0, </?) and $; (f, 7r, </j) that is 
on x 3 axis. 

The angle spherical bivector functions is defined as ©j m ($, ip) = W- m §™ 
[7], where ^ m = Z*- m (cosi?) are zonal functions, §™ = §™(i9, y?) are secto- 
rial functions. According to [12] we have for the zonal functions 



i} m (i) = c , s} m (-i) = e m * 2 ' = e m (-i)' , (7.5) 

where I is considered integer. According to (8.8 in [7]) §™(0, </?) = CD™ - , 
g«(7r, = -(K™. for j = 0, ±1. As result we have 

©f (0, y?) = $™ GDf , ©; m (7r, </?) = (-l) m G CD!! 1 , , (7.6) 
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where j = 0, ±1. 

I consider both with progressing plane wave as with standing 

one. Thus I take expression (16.51) or (I6.17P for Y accordingly. For Y I 
take representation (15. ip with coefficients (17. 2p . where a number (^ l ± m {u) q ) 
is calculated with the help of relations (I7.1bl) with (15.261) and (16.131) or 



( I6.20p . I substitute these expressions for Y, Y into (17.4bl) . (I7.4cp and use 
( 17.61) . (16.6bp . As result we have for the progressing wave 



{€} =(|© + | 2 + |©_| 2 )/4vr , (7.7a) 

$ l £(r,0,<p) = ( 2l i + l)^- 1 (±\&-\ 2 - \&+\ 2 ) tan(^ ± ) , (7.7b) 

<f> l l(r,ir,ip) = (7.7c) 

and for the standing wave 



[6) = {\& + \ 2 + m 2 ) /2n , (7.8a) 

* l }(r,0,<P) =^(r,7r,^) = (2/, + r)^- 1 (|ffi_| 2 - |© + | 2 ) tan(^) , 

(7.8b) 

$ l l(r,0,<p) =$ l l(r,*,tp) = -(2l q + l)Cu- 1 *- 1 [(|©_| 2 + \& + \ 2 ) 

+ {©+ &_ + *&+ *©_)] tan(g^) , (7.8c) 

where we must take tan(^^) ^ or tan(€?i) ^ 0. 

We can see in formulas (17. 7p and (17. 8p the following. The function 
a; 3 axis is asymmetrical for the case with progressing wave 
( 17.71) and it is symmetrical for the case with standing wave (17. 81) . The val- 
ues of the function $+(f, cp) at a; 3 for the standing wave can be obtained 
by summation of the appropriate values for progressing waves with oppo- 
site directions of propagation (compare (]7.7bj) and (j7.7cj) with (I7.8bl) ). But 
the values of the function (f, $, <p) at a; 3 can not be obtained by the sum- 
mation (compare (I7.7bp and (I7.7cl) with (I7.8c[) ) . The function <&±(f, (p) 



can take positive or negative values. It should be noted that here we have 
the behavior which is typical for interference of waves (Y and Y for this 
case) . 

The gradient of the averaged energy density creates an appropriate co- 
directed force for another electromagnetic particles and distorts rays of 
light. This effect appear because the effective metric component contains 
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the averaged energy density [21 El E] : 

moo = m 00 - 4vr x 2 {£} 



(7.9) 



As we see in (17.71) and (17. 8p . the near resonance wave solution with 
one acting plane wave can give both attraction ($± > 0) an d repulsion 
< 0) effects. Also there is an angular dependence. But it would 
appear reasonable that the angular dependence can be canceled with taking 
into consideration of all acting plane waves contained in the equilibrium 
wave background, as well as by an averaging for a bulk of electromagnetic 
particles just as in macroscopic bodies. The appropriate averaged energy 
density can have the form (I7.4ap with substitution $± — > $ = Constant. 
In principle, the following three cases is possible: $ > 0, $ = 0, and 
$ < 0. The case $ > corresponds to attraction which is essential for 
explanation of gravitation. However, in the context of such approach we 
must not exclude a possibility of special and artificial field configurations 
which give the repulsion or antigravitation effect, for example, with some 
coherent states of matter. 

8 Wave solution at zero 
and beyond linearization 

The radial functions C+ w and ^ contained in solution (15 . 1 2[) have different 
behavior near zero. The functions £+ w has the main term proportional 
to f -2 (I5.15ap . just like basic solution (14. ip . But the main term of the 
functions is proportional to f ( 15.15bl) . 

If we have the resonance increase for the amplitudes @ l ™ of the solution 
(see ( I6.13P and (16.201) ) then linearization (13.41) is no longer well approxi- 
mation. But the plus component resonance & l + will rather destroy the 
linearization then the minus component one because of the different 
behaviour for the components near zero. 

Thus we must regard the minus components resonance frequencies uo~ 
to be more plausible then the plus ones (see Tab. |2]). 

In any case, to avoid the infinite amplitude @ l ± at a resonance frequency 
we must consider a nonlinear resonance by taking into account highest 
terms in expansion (13.41) or solving the problem with original constitutive 
relation (121I . 
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It is known, in general resonance phenomena theory for oscillating sys- 
tem, that resonance frequencies for the case of nonlinear resonance most 
often depend on an influence amplitude. But this dependence can be weak 
or strong. What is more, we can suppose that some resonance influences 
will destroy the singularity under consideration and will give birth for sin- 
gularities of another types. 

It should be noted also that a veritable wave solution of the suitable 
nonlinear problem can essentially differ from the wave solution of the lin- 
earized problem. 

9 Conclusions 

The present investigation gives the method and solutions concerning the 
problem linearized around the static dyon solution for Born-Infeld non- 
linear electrodynamics. In particular, the wave resonance effect on static 
electromagnetic particle is discovered. This effect appears explicitly, if the 
the static electromagnetic particle is imbedded to an equilibrium electro- 
magnetic wave background. 

But we must consider an electromagnetic wave background to be always 
in space. Thus we can not have a static solution without the appropriate 
wave modes. We can suppose that the wave modes will essentially influence 
on behavior and manifestations of the particle solution in its interaction 
with another particles. The demonstrated possible connection with the 
gravitational interaction is one of such manifestations. 
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A Main designations 



Symbol 


Name 


Appearance 




Metric tensor, effective metric 


(123)1. 07^]) 


<S>, A 


^^rin TYi^"i"Ti r*£i 1 q ti H a cvtti tti f^Ti r*fi 1 t»t*oH i l r*"i" c 
kj y iiiiiic vl lKjCXiL cliiu. do y iiiiiic li ii^cii ui uulil-id 


AO 1 h 




Basis vectors and bivectors 


m 


E, Ei, H, Hi 


Intensities of fields: electric and magnetic 


(I2.20 


D, D\ B, 


Inductions: electric and magnetic 


02. 2p 


Y, Z 


Quasi-bivectors of electromagnetic field 


(12.20 


*C 


Hyperconj ugation 


02. 41) 




Energy density and momentum density 


(1X7]) 


Y, Y, Y 


Basic solution for linearization and lin- 
earized solution 


([3~T]). (|63]) 


p C e g 


Electromagnetic, normed, electric, and 
magnetic charges 


(14 11) 


f, f, cD 


Radius of SDS, dimensionless radius and 


(I4.ip. 05.111) 




frequency 




Angle spherical functions: general, secto- 
rial, and zonal 


05. ip. (I7.5P 




OUlloLcLlllo 1UI tilt: llUcdl 1/icU oUlULlUIl 


(lo.lzll. do. loll 


S± 


Radial functions of the problem 


(|5.12p 


O r_1 o r_1 


Infinitesimals of the same and higher order 


05.151). 07.21) 




TT 111 * i 

Hyperreal and hypenmagmary parts 




05.171). 05.281) 




Radial ^nVipriral fimptinns' 

1 111 VI 1(11 1 711V 1 1 V < 11 lUlUjllVvllOi 


i fih 9sn 

IIU.1UII. IIO.ZOII 


|C|, (<7> 


Modulus and argument of hyperscalar 


05.221) 




Angle cylindrical functions 


06. 5p 




Polarization coefficients 


06. 5p. 06.171) 




Spherical coefficients of plane wave 


06.6b) 




Kronecker symbol 


(1631) 




Response functions 


16.13b). 07.1b! 








Averaging operation 
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